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Abstract. Radial motions of stars in nearby moving clusters are determined from accurate proper motions 
and trigonometric parallaxes, without any use of spectroscopy. Assuming that cluster members share the same 
velocity vector (apart from a random dispersion), we apply a maximum-likelihood method on astrometric data 
from Hipparcos to compute radial and space velocities (and their dispersions) in the Ursa Major, Hyades, Coma 
Berenices, Pleiades, and Praesepe clusters, and for the Scorpius-Centaurus, a Persei, and 'HIP 98321' associations. 
The radial motion of the Hyades cluster is determined to within 0.47 km s -1 (standard error), and that of 
its individual stars to within 0.6 km s -1 . For other clusters, Hipparcos data yield astrometric radial velocities 
with typical accuracies of a few km s _1 . A comparison of these astrometric values with spectroscopic radial 
velocities in the literature shows a good general agreement and, in the case of the best-determined Hyades cluster, 
also permits searches for subtle astrophysical differences, such as evidence for enhanced convective blueshifts 
of F-dwarf spectra, and decreased gravitational redshifts in giants. Similar comparisons for the Scorpius OB2 
complex indicate some expansion of its associations, albeit slower than expected from their ages. As a by-product 
from the radial-velocity solutions, kinematically improved parallaxes for individual stars are obtained, enabling 
Hertzsprung-Russell diagrams with unprecedented accuracy in luminosity. For the Hyades (parallax accuracy 
0.3 mas), its main sequence resembles a thin line, possibly with wiggles in it. Although this main sequence has 
underpopulated regions at certain colours (previously suggested to be 'Bohm-Vitense gaps'), such are not visible 
for other clusters, and are probably spurious. Future space astrometry missions carry a great potential for absolute 
radial-velocity determinations, insensitive to the complexities of stellar spectra. 
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1. Introduction 

This paper is the third of a series on the determina- 
tion of stellar radial motion by purely geometric measure- 
ments. Such astrometric radial velocities allow to disen- 
tangle stellar motion from other astrophysical phenomena 
causing spectroscopic line shifts, such as internal motions 
in stellar atmospheres and gravitational redshifts. Paper I 
(Dravins et al. 1999b| ) discussed different methods which 
permit radial velocities to be determined independent of 
spectroscopy. Among these, the moving-cluster method of 
changing angular separation permits radial-velocity accu- 
racies on a sub-km s~ x level to be reached already with 
current astrometric measurements, such as those from the 



Hipparcos space mission (ESA 1997). That method anal 
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yses the changing angular extent of a star cluster as it 
approaches or recedes from the Sun, assuming that the 
member stars share the same average velocity vector rela- 
tive to the Sun. The radial-velocity component makes the 
cluster appear to expand or contract due to its changing 
distance. This relative rate of apparent contraction equals 
the relative rate of change in distance to the cluster, from 
which a linear velocity (in km s" 1 ) follows whenever the 
absolute distance is known from trigonometric parallaxes. 



That it is possible, in principle, to determine radial 
velocities from astrometry has been known for a long time. 
Attempts were made e.g. by Petrie (194S, 1963| ) and Eggen 
(1984, 1998) to derive radial velocities that are at least 
partially independent of spectroscopy. The availability of 
the Hipparcos results allowed to investigate the kinematics 
and memberships of moving groups and clusters in great 
detail, which in turn made it possible to apply our moving- 
cluster method to reach better accuracies. The detailed 
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Fig. 1. Map of the full sky, showing those stars in clusters and associations, whose radial velocities were astrometrically 
determined from Hipparcos data. Symbol shape identifies different clusters; symbol size denotes apparent magnitude 
Hp (~ my), while symbol shading denotes B — V (note how some clusters are dominated by very blue stars). The 
Aitoff projection in equatorial coordinates is used, with S — 0° on the major axis and a = 180° on the minor axis. 
Right ascension (a) increases to the left. 



mathematical formulation of the method was presented in 



Paper II (Lindegren et al. 2000) 



As a by-product of the moving-cluster method, the 
distance estimates to individual cluster stars are often 
significantly improved compared with the original [here 
Hipparcos] parallax measurements. These kinematically 
improved parallaxes (Paper II) result from a combina- 
tion of trigonometric and kinematic distance information, 
where kinematic distances follow from the observed proper 
motions and the derived cluster velocity. The improved 
distances sometimes allow the Hertzsprung-Russell dia- 
gram of a cluster to be studied with unprecedented reso- 
lution in absolute magnitude. 

In this paper we apply this 'moving-cluster method' to 
several nearby open clusters and associations. Astrometric 
radial velocities and kinematically improved parallaxes are 
deduced for stars in the the Ursa Major, Hyades, Coma 
Berenices, Pleiades, and Praesepe clusters, and for the 
Lower Centaurus Crux, Upper Centaurus Lupus, Upper 
Scorpius, a Persei, and 'HIP 98321' associations. 

2. Potential of the moving-cluster method 

The achievable accuracies of the moving-cluster method 
were discussed in Paper I. We recall that the best radial- 
velocity accuracy is obtained for star-rich nearby clusters 
with large angular extent, large proper motions, small in- 



ternal velocity dispersions, and small rates of cluster ex- 
pansion. The improved astrometric accuracies expected 
from future space missions will somewhat lessen these con- 
straints, although the intrinsic limitations set by internal 
velocities cannot be overcome by increasing observational 
accuracy. 

Only about ten clusters and associations can be mean- 
ingfully studied already with current astrometric accura- 
cies in the milliarcsecond range. Furthermore, good astro- 
metric data are available only for their relatively brighter 
stars. The distribution on the sky and the geometries of 
these nearby clusters are shown in Figs, [l] and |[ together 
with their stellar populations, and their proper-motion 
patterns. Only a few among these (Hyades, Pleiades, 
Coma Berenices, Praesepe) have a utilizable stellar pop- 
ulation spanning many spectral types; several of the oth- 
ers are heavily dominated by early-type stars. The main 
reason for this is of course the limiting magnitude of the 
Hipparcos mission. The areas subtended on the sky differ 
greatly: some OB-associations spread out over much of a 
hemisphere, while some other clusters are very localised. 
Although the great spatial extents of the Ursa Major clus- 
ter and the Scorpius-Centaurus associations in principle 
are advantageous for the accuracy obtainable, the partly 
unknown expansion rates and internal velocity patterns 
of these younger stellar groups actually limit the accuracy 
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Fig. 2. Proper motions of the programme stars over 200,000 years. Best radial-velocity accuracy is obtained in rich 
nearby clusters with large angular extent, and large proper motions. However, the accuracy in the largest associations 
(Ursa Major, Scorpius-Centaurus) is limited by the partly unknown expansion of these systems. Stellar paths in the 
Ursa Major group (shown shaded) cover large areas of the sky. The thickness of the proper-motion vectors is inversely 
proportional to stellar distance: the closest star is Sirius and the two next ones are faint red dwarfs. Proper motions 
vary greatly among different clusters. 



in their radial-velocity determinations. For more data on 
these clusters, see Table 4 in Paper I. 

3. Exploitation of the moving-cluster method 

3.1. Basic cluster model 

The mathematical procedure described in Paper II yields 
maximum-likelihood estimates for the space velocity of the 
cluster centroid vq, for the internal velocity dispersion of 
stars within the cluster d v , and for the kincmatically im- 
proved parallax of each star 5?,, i = 1, 2, ... n. (We reserve 
index for the centroid. The caret ~ signifies an estimated 
value.) Although additional model parameters could be in- 
cluded, e.g. to describe a possible rotation or non-isotropic 
dilation of the cluster, the present studies are restricted 
to the 'basic cluster model' in which no such systematic 
velocity patterns are assumed. In Paper II it was shown, 
through Monte Carlo simulations of the Hipparcos obser- 
vations of the Hyades cluster, that the presence of any 
reasonable amount of rotation and shear in the actual 
cluster will not significantly bias the solution for the cen- 
troid velocity, even though the analysis is restricted to the 
basic model. One important exception concerns the possi- 
ble isotropic expansions of gravitationally unbound asso- 
ciations. These could indeed introduce significant biases, 



3.2. Observational data 



which are discussed separately in Sect. 5.7 



The Hipparcos Catalogue (ESA |l997 ) provided input data 
for each star i in the form of positions in barycentric right 
ascension OLi, declination Si, trigonometric parallax 7?,, the 
proper- motion components Jl a i and Jlgi , and standard de- 
viations and correlation coefficients for the latter three. 
(The tilde ~ signifies an observed value, the uncertainty 
of which needs to be taken into account in the estima- 
tion procedure.) The positions and proper motions in the 
Hipparcos Catalogue are referred to the barycentric ICRS 
reference system, and consequently all resulting velocities 
are also in that system. 

For each cluster or association, an initial sample of 
probable member stars was identified from the literature, 
mainly from studies based on Hipparcos data. However, 
since the mathematical formalism for obtaining radial ve- 
locities is strictly applicable only to stars sharing the 
same average velocity vector (with a random spread about 
that value), cluster non-members and binary stars in non- 
modelled orbits must be removed from the sample as 
far as possible. This was done using an iterative rejec- 
tion procedure described in Paper II. Monte Carlo simu- 
lations (Sect. 4.2 in that paper) showed that this proce- 
dure works best when adopting the goodness-of-fit rejec- 
tion limit gn m — 15. As illustrated there for the Hyades, 
this gave the lowest scatter in the centroid radial veloc- 
ity, as well as in other quality indicators. Therefore, unless 
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Table 1. Estimated space- velocity components and internal velocity dispersions of clusters and associations analyzed 
with the moving-cluster method, using astrometric data from the Hipparcos main catalogue. n acc is the number of 
stars retained in the cleaned sample for each cluster, n re j is the number of stars removed in the cleaning process; Vq x , 
Vq v , and v 0z are the equatorial (ICRS) components of the estimated space velocity of the cluster centroid; a v is the 
estimated internal velocity dispersion among individual stars, calculated as described in Paper II, Appendix A. 4. The 
last three columns give the equatorial coordinates (cto, So) [deg] of the adopted centroid of the cluster (or, in the case of 
Ursa Major, its core) and the space velocity component «or toward that direction, i.e. an approximate radial velocity 
of the cluster as a whole. Uncertainties are given as ±1 standard error. All velocities are in km s _1 . The electronic 
version of the table contains, additionally, the six equatorial components of the formal covariance matrix Cov(vo), the 
spherical equatorial coordinates of the convergence point with standard errors, and the total velocity with its standard 
error. 



Name 


^acc 




Vox 




v 0y 




vo 


Z 












So 






Ursa Major 


77 


4 


+8.44 


± 


0.41 


-12.19 


± 


0.39 


-10.16 


± 


0.43 


2.82 


± 


0.23 


187.3 


+56.4 


-12.24 


± 


0.46 


Hyades 


168 


29 


-5.90 


± 


0.13 


+45.65 


± 


0.34 


+5.56 


± 


0.10 


0.49 


± 


0.04 


66.5 


+16.9 


+39.42 


± 


0.36 


Coma Berenices 


40 





-0.82 


± 


0.96 


+4.57 


± 


0.15 


-4.11 


± 


0.48 


0.47 


± 


0.09 


187.5 


+26.4 


-1.64 


± 


1.07 


Pleiades 


60 





+1.99 


± 


2.20 


+22.95 


± 


3.34 


-18.73 


± 


1.82 


0.50 


± 


0.13 


56.4 


+24.0 


+10.85 


± 


4.36 


Praesepe 


24 





-1.46 


± 


9.03 


+48.1 


± 


10.9 


+2.00 


± 


5.02 


0.67 


± 


0.23 


130.2 


+19.6 


+36.2 


± 


15.0 


Lower Cen Crux 


179 


1 


-0.94 


± 


0.29 


+18.36 


± 


0.15 


-8.59 


± 


0.46 


1.13 


± 


0.07 


189.6 


-56.2 


+5.95 


± 


0.53 


Upper Cen Lupus 


218 


3 


-4.01 


± 


0.27 


+16.47 


± 


0.32 


-12.91 


± 


0.37 


1.23 


± 


0.08 


230.3 


-41.6 


+1.01 


± 


0.51 


Upper Scorpius 


120 





-3.73 


± 


0.56 


+9.36 


± 


1.09 


-14.57 


± 


0.57 


1.33 


± 


0.12 


243.4 


-24.1 


-0.17 


± 


1.33 


Sco OB2 


510 


11 


-1.72 


± 


0.15 


+18.19 


± 


0.15 


-10.43 


± 


0.21 


1.52 


± 


0.06 


225.1 


-43.9 


-1.17 


± 


0.26 


a Per (Per OB3) 


78 


1 


-3.23 


± 


0.89 


+27.15 


± 


1.26 


-11.76 


± 


1.64 


0.71 


± 


0.13 


52.9 


+47.8 


+4.53 


± 


2.18 


'HIP 98321' 


59 





-3.45 


± 


0.66 


+15.55 


± 


1.26 


-12.27 


± 


1.13 


2.56 


± 


0.26 


297.5 


+39.4 


-19.68 


± 


1.74 



otherwise stated, all samples discussed here were cleaned 
according to this criterion. 

Solutions were primarily obtained using data from the 
main Hipparcos catalogue but, for some clusters, w e use d 
also data from the Tycho-2 catalogue (H0g et al. 20(Xf ). 
The accuracies in the latter are generally somewhat worse, 
but since its proper motions incorporate about a cen- 
tury of ground-based observations, the segregation of long- 
period binaries may be improved. 

3.3. Clusters and associations studied 

Table 4 in Paper I lists some 15 clusters and associa- 
tions, for which Hipparcos-type accuracies could poten- 
tially yield viable astrometric radial velocities, i.e. with 
standard errors less than a few km s . In practice the 
resulting accuracies depend on several additional factors, 
not considered in that survey, such as the number of mem- 
ber stars actually observed by Hipparcos, the position of 
the cluster on the sky, the statistical correlations among 
the astrometric data, and the procedures used to obtain 
a clean sample. The simplest way to find out whether our 
method 'works' on a particular group of stars is to make 
a trial solution. We have done that for all the potentially 
interesting clusters and associations, based on the list in 
Paper I supplemented with data from the compilations 
in the Hipparcos Input Catalogue (Turon et al. 1992), by 



Robichon et al. (|999|), and by de Zeeuw et al. ( |1999| ). 

The original criterion for including a cluster or asso- 
ciation in the present study was that it yielded a valid 
solution with the basic cluster model (Sect. 3.1), includ- 



turned out to be the case for four clusters (Ursa Major, 
Hyades, Coma Berenices, and Praesepe) and four associa- 
tions (Lower Centaurus Crux, Upper Centaurus Lupus, 
Upper Scorpius, and 'HIP 98321'). For two more, the 
Pleiades cluster and the a Persei association, reasonable 
solutions could be obtained by assuming zero velocity dis- 
persion in the maximum-likelihood procedure (note that 
the dispersion could still be estimated from the proper- 



motion residuals, as explained in Sect. 3.4). Considering 
the astrophysical importance of these clusters, they were 
therefore included in the study. A separate solution was 



also made for the Sco OB2 complex (Sect. 5.4) 



ing a non-zero estimate of the velocity dispersion. This 



3.4. Mathematical bias and noise in the solutions 

Resulting space velocities and internal velocity dispersions 
for the 11 clusters and associations are in Table |l]. vq x , 
Voy, and vq z are the equatorial components (ICRS co- 
ordinates) of the estimated space velocity of the cluster 
centroid, while a v is the velocity dispersion in each co- 
ordinate, i.e. the standard deviation of peculiar velocities 
along a single axis. 

The maximum-likelihood estimation tends to under- 
estimate the velocity dispersion, as examined through 
Monte Carlo simulations in Paper II. In Appendix A. 4 of 
that paper we gave an alternative procedure to estimate 
the velocity dispersion from the proper-motion residuals 
perpendicular to the centroid velocity projected on the 
sky. This was shown to give nearly unbiased results. The 
velocity dispersions d v in Table ^ and elsewhere in this pa- 
per have therefore been estimated through this alternative 
procedure. 



S0ren Madsen et al.: Astrometric radial velocities III. 



5 



As was also described in Paper II, the radial-velocity 
errors among individual stars in the same cluster are not 
statistically independent, but may carry a significant pos- 
itive correlation. For each star, the error contains a [nearly 
constant] component being the uncertainty in the cluster 
velocity as a whole, plus a random component correspond- 
ing to the physical velocity dispersion among the individ- 
ual stars. Averaging over many stars in a given cluster 
averages away the influence of the velocity dispersion, but 
has only little effect on the error in the radial velocity 
of the cluster centroid. This quantity, discussed already 
in Paper I (e.g. its Table 4), is therefore a limiting ac- 
curacy in the average astrometric radial velocity of stars 
in any one cluster. For certain applications, effects of this 
noise can be lessened by averaging over different clusters 
(whose errors are not correlated), e.g., when searching for 
systematic differences between astrometric and spectro- 
scopic radial-velocity values. 

It should be noted that the selection process used to 
arrive at the final sample may have a significant impact 
on the estimated internal velocity dispersion. The cleaning 
process successively removes those stars that deviate most 
from the mean cluster velocity, thus successively reducing 
a v for the 'cleaner' samples. While designed to remove 
non-members and other outliers, this procedure naturally 
affects also the mean characteristics of the remaining stars. 
For example, in the case of the Hyades, stars in the out- 
skirts of the cluster are preferentially removed during the 
rejection procedure, meaning that the resulting clean sam- 
ples more or less correspond to the stars within the tidal 
radius. In the case of the OB associations, we obtain ve- 
locity dispersions of about 1 km s _1 , a factor of two lower 
than th e esti mates for th e Orio n Nebula Cluster (Jones & 
Walker 1988j ; Tianet al. 1996 ) and other nearby associa- 
tions (Mathieu 1986 ). Thus, although we believe that the 
velocity dispersions reported here correctly characterize 
the retained samples, they are not necessarily representa- 
tive for the cluster or association as a whole. 

3.5. Calculation of astrometric radial velocities 

3.5.1. The stringent definition of 'radial velocity' 

Recognizing the potential of astrometric radial velocities 
determined without spectroscopy, a resolution for their 
stringent definition was adopted at the General Assembly 
of the International Astronomical Union held in 2000. This 



resolution (Rickman 2001 ) defines the geometric concept 
of radial velocity as v r = db/dtB, where b is the barycen- 
tric coordinate distance to the object and the barycen- 
tric coordinate time (TCB) for light arrival at the solar 
system barycentre. This definition is analogous to the con- 
ventional understanding of proper motion as the rate of 
change in barycentric direction with respect to the time 
of light reception at the solar-system barycentre. 

In this work, we follow this IAU definition of 'astro- 
metric radial velocity'. The difference with respect to al- 
ternative possible definitions is on the order of vf. / c, with 



c = speed of light (Lindegren et al. 1999 Lindegren & 
Dravins, in preparation). Most population I objects (in- 
cluding all clusters and associations considered in this pa- 
per) have low velocities, \v r \ < 50 km s _1 , resulting in 
only very small differences, < 10 m s _1 , between possible 
alternative definitions. 



3.5.2. Radial velocities for individual stars 

In the basic cluster model, the estimated radial velocity 
of an individual star is given by 



r'iV Q , 



(1) 



where is the unit vector towards star i and v$ is the 
estimated space velocity of the cluster as a whole (actually 
of its centroid). The prime (') denotes the transpose of the 
vector. In terms of the equatorial coordinates (oii,8i) of 
the star we have 



Vri — vq x cos Si cos Oi + vo y cos Si sin o>i + vq z sin Si 



(2) 



where (vq x , vq v , Vq z ) are the equatorial velocity compo- 
nents as listed in Table |l| (or as determined by other 
means). We emphasize that Eq. (FJ) applies to any star 
that shares the cluster motion, irrespective of whether 
that star was present in the database used to determine 
the cluster motion in the first place. 

The standard error e(v r i) of the individual radial ve- 
locity is computed from 



e(«ri) 2 = r ■Cov(v )r i + '. 



(3) 



where Cov('Uo) is the 3x3 submatrix in Cov(0) of all 
model parameters, referring to the centroid velocity [cf. 
Eq. (A18) in Paper II]. The first term in Eq. (FJ) represents 
the uncertainty in the radial component of the common 
cluster motion, while the second represents the contribu- 
tion due to the star's peculiar motion. 

The complete covariance matrix Cov(vo) is only given 
in the electronic (extended) version of Table hi The printed 
Table |l| gives (following the ± symbol) the standard er- 
rors e(vox) etc. of the vector components; these equal the 
square roots of the diagonal elements in Gov(vq). Also 
given for each cluster is the standard error of the radial 
component vq t — r' VQ of the centroid motion. This was 
computed from 



e(vor) 2 = r' Cov(v )r , 



(4) 



where rQ is the unit vector towards the adopted centroid 
position (tto; 8q) specified in the table. vo r can be regarded 
as an average radial velocity for the cluster as a whole, and 
its standard error (squared) can be regarded as a typical 
value for the first term in Eq. (|). Thus, for any star not 
too far from the cluster centroid, the total standard er- 
ror of its astrometric radial velocity can be approximately 
computed as [e(vo r ) 2 + c? 2 ] 1 / 2 , using only quantities from 
the printed Table |l|. 
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3.6. Kinematically improved parallaxes 

Our maximum-likelihood estimation of the cluster space 
motions also produces estimates of the distances to all 
individual member stars. A by-product of this moving- 
cluster method is therefore that individual stellar dis- 
tances are improved, sometimes considerably, compared 
with the original trigonometric determinations. This im- 
provement results from a combination of the trigonomet- 
ric parallax 7r tr ig with the kinematic (secular) parallax 
TTkin = An/v t derived from the star's proper motion fj, 
(scaled with the astronomical unit A) and tangential ve- 
locity vt , the latter obtained from the estimated space ve- 
locity vector of the cluster. The calculation of secular par- 
allaxes and kinematic distances to stars in moving clus- 
ters is of course a classical procedure; what makes our 
'kinematically improved parallaxes' different from previ- 
ous methods is that the values are derived without any 
recourse to spectroscopic data (for details, see Papers I 
and II). 



Dc Bruijnc ( 1999b ) applied the present method to the 
Scorpius OB2 complex in order to study its HR diagram 
by means of the improved distances. His 'secular paral- 
laxes' are essentially the same as our 'kinematically im- 
proved parallaxes', being based on the same original for- 
mulation by Dravins et al. (1997). The main differences 
are in the choice of rejection criteria (de Bruijne uses 
5iim = 9 versus our 15) and in the practical implementa- 
tion of the solution (downhill simplex versus our use of an- 
alytic derivatives). De Bruijne also made extensive Monte 
Carlo simulations which demonstrated that the distance 
estimates are robust against all systematic effects consid- 
ered, including cluster expansion. For the accuracy of the 
secular parallaxes, de Bruijne ( 1999b| ) used a first-order 
formula (his Eqs. 16 and 17) which explicitly includes a 
contribution from the (assumed) internal velocity disper- 
sion, but neglects the contribution from the trigonometric 
parallax error (cf. Eq. 11 in Paper I). By contrast, our 
error estimates are derived directly from the maximum- 
likelihood solution (Paper II, Appendix A. 3), which in 
principle takes into account all modelled error sources but 
in practice underestimates the total error as discussed be- 
low. As a result, our error estimates are somewhat smaller 



than those given by de Bruijne ( 1999b ). 

The standard errors for the estimated parallaxes given 
in this paper are the nominal ones obtained from the 
maximum-likelihood estimation, which could be an un- 
derestimation of the actual errors. Determination of real- 
istic error estimates would require extensive Monte-Carlo 
experiments based on a detailed knowledge of the actual 
configuration of stars, their kinematic distributions, etc. 
This information is in practice unavailable except in ide- 
alised simulations, and we therefore choose not to intro- 
duce any ad hoc corrections for this. As an example of 
the possible magnitude of the effect, the Hyades simula- 
tions in Paper II could be mentioned: in that particular 
case, the nominal standard errors required a correction by 



a factor 1.25 to 1.28 in order to agree with the standard 
deviations in the actual sample. 

4. Radial velocities for stars in open clusters 

Partial data for astrometric radial velocities and kine- 
matically improved parallaxes of individual stars in more 
nearby clusters, obtained from Hipparcos data, are in 
Table |2[ However, the complete listing for the more than 
1000 stars in all clusters and associations (some with also 
Tycho-2 solutions), is only in the electronic Table |[ 

For some clusters, Fig. || shows a comparison between 
spectroscopic radial velocities compiled from the litera- 
ture and our currently determined astrometric values. In 
all cases (also the ones not shown here), the astrometric 
values agree with the spectroscopic ones within the er- 
ror limits, verifying the consistency of our method. For 
several stars (especially rapidly rotating early-type ones 
with smeared-out spectral lines), the astrometric accura- 
cies are actually substantially better than what has been 
possible to obtain from spectroscopy. The high accuracies 
realize d for the Hyades enable more detailed comparisons 
(Sect. |4.2.l|) . 

A subset of known or suspected binary stars for which 
Hipparcos measurements could be perturbed are defined 
as stars that are either visual binaries with magnitude dif- 
ference Am < 4 mag and a separation p < 20 arcsec ac- 
cording to HIC, the Hipparcos Input Catalogue (Turon et 
al. 1992 ); known spectroscopic binaries; or flagged as sus- 
pected ones (identified in the Hipparcos Catalogue as a so- 
lution of type component, acceleration, orbital, variability- 
induced mover, or stochastic). In Fig. || and later, these 
stars are plotted in gray. 

For finite astrometric accuracy, distant clusters of 
small angular extent basically yield a single astrometric 
velocity for all stars. This effect is seen in Fig. || for Coma 
Berenices and the Pleiades, and is similar also for Praesepe 
and the 'HIP98321' clusters (not plotted). 

4.1. The Ursa Major cluster 

The initial sample consisted of 81 stars, being th e sum 
of the compila tions from Soderblo m fc Mayor ( 1993 ), 
Dravins et al. ( ]1997| ) and Montes ( |2000D . A few of the 
stars identified in a previous kinematic search (Dravins 
et al. 1997 ) lacked spectroscopic radial velocities from 
the literature. Spectroscopic observations of these stars 
by Gullberg & Dravins (private comm.) made with the 
ELODIE spectrometer at Observatoire de Haute-Provence 
confirmed three of them as probable members. The rejec- 
tion procedure removed four stars, producing a final sam- 
ple of 77. 

We chose to include stars not only from the core but 
also from the extended halo (the moving group sometimes 
called the Sirius stream), in order to improve the statisti- 
cal weight of the solution. Not surprisingly, this led to a 
relatively high velocity dispersion among individual stars. 
Ursa Major may thus be viewed as a dissolved cluster 
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Fig. 3. Astrometrically determined radial velocities com- 
pared with spectroscopic values from the literature, for 
stars in four open clusters. The diagonal lines follow the 
expected relation iv(astrom) ~ iv(spectr). Black symbols 
denote single stars while certain or suspected binaries are 
in grey. The top three frames are kinematic solutions ob- 
tained from Hipparcos data only. In some cases, includ- 
ing the Pleiades (bottom) , somewhat better accuracies are 
reached using data from the Tycho-2 catalogue, which in- 
corporates almost a century of proper-motion data. 



moving under influence of the Galactic gravity field. Such 
a fate may be normal for looser open clusters reaching the 
age (300 Myr) of Ursa Major (Soderblom & Mayor |1993| ). 
For convenience we use the designation 'cluster' for the 
entire sample of stars. 

One could imagine that moving groups like Ursa Major 
could be ideal targets for the method since they have 
great angular extents. Unfortunately, as Ursa Major il- 
lustrates, the high velocity dispersion causes the errors 
of the estimated astrometric radial velocities to be large. 
However, it should be noted that the core stars have 
a much lower velocity dispersion. For instance, we get 
a v — 1.05 ±0.22 km s _1 for the 13 co re (nucleu s) sta rs de- 
fined by Soderblom & Mayor fll993| ). Wielen ([L978D con- 
sidered six core stars with very well determined proper 
motions and found the velocity dispersion to be on the or- 
der of 0.1 km s _1 (we get 0.09 ± 0.03 km s _1 for the same 
six stars), adding that the much larger stream, or moving 
group, has a dispersion of ~ 3 km s , in good agreement 
with our value. Since we assume this larger dispersion for 
all the stars, the standard errors are probably overesti- 
mated for the core stars. 

In the case of Ursa Major, no real gain results from the 
kincmatically improved parallaxes, simply because the rel- 
ative accuracy in the Hipparcos parallax values is already 
very good, given the proximity of this cluster. 

The spectroscopic velocity values used for Fig. |^ were 
taken from Soderblom & Mayor (1993), Duflot et al. 
(1995), and in a few cases, also the ELODIE observations 
by Gullberg & Dravins (private comm.). In the case of 
Duflot et al., no explicit error is quoted, but rather a flag 
which seems to correspond to a numerical value that is 
used in the Hipparcos Input Catalogue: those values were 
adopted here. 



4.2. The Hyades cluster 

The Hyades cluster (Melotte 25) is the classic example 
of a moving cluster. Its kinematic distance, derived from 
a combination of proper motions and spectroscopic radial 
velocities, has been one of the fundamental starting points 
for the calibration of the photometric distance scale (e.g. 
Hanson |1975| , Gunn et al. |1988| , Schwan |l99lj and ref- 
erences therein). Of course, the recent availability of ac- 
curate trigonometric parallaxes has now superseded this 
method for distance determination. 

The first detailed study of the distance, structure, 
membership, dynamics and age of the Hyade s clus ter, us- 
ing Hipparcos data, was by Perryman et al. ( 1998 ). From 
a combination of astrometric and spectroscopic radial ve- 
locity data, using 180 stars within a radius of 20 pc, 
they derived the space velocity for the cluster centroid, 
v = (-6.32, +45.24, +5.30) km s -1 . They also estimated 
that the true internal velocity dispersion, near the centre 
of the cluster, is in the range 0.2 to 0.3 km s _1 . 

Our initial sample of stars was selected from the fi- 



nal membership assigned by Perryman et al. (1998), viz. 
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197 stars classified as probable members (S = 1 in their 
Table 2): this equals our sample HyO in Paper II. 

In earlier kinematic studies of the Hyades, systematic 
errors in the proper motions have been of major concern, 
and the probable cause of discrepant distance estimates. In 
our application, the solution is also sensitive to such errors, 
but we expect that the high internal consistency of the 
Hipparcos proper motion system, and its accurate linking 
to the (incrtial) extragalactic reference frame (Kovalevsky 
et al. 1997), have effectively eliminated that problem. 

Using the methods described in Paper II, the astro- 
metric radial velocity and its standard deviation for each 
individual star is obtained through Eqs. (||) and (||), even 
for stars not retained in the final sample. Individual par- 
allaxes follow directly from the estimation procedure, see 
Table § 

The analysis was repeated using proper motions from 
Tycho-2 instead of the Hipparcos Catalogue. Being based 
on observations covering a much longer time span, the 
Tycho-2 data are expected to be more precise for binaries 
with periods from a few years to ~ 100 years. By co- 
incidence, the procedure of cleaning the sample rejected 
the same number of stars as in the Hipparcos case, al- 
beit different ones (electronic Table |2|) . Compared to the 
Hipparcos solution (0.49 ± 0.04 km s _1 ), we get a smaller 
velocity dispersion using Tycho-2 (0.34 ± 0.03 km s _1 ). It 
can be noted that Makarov et al. ( 2000] ) investigated this 
dispersion in the Hyades as a test of the Tycho-2 proper 
motions, finding a dispersion close to our latter value. The 
two solutions for the cluster centroid velocity are equal to 
within their uncertainties. However, in the radial direc- 
tion there is a systematic difference of ~ 0.9 km s _1 , 
in the sense that the astrometric radial-velocity values 
from Tycho-2 are smaller than those from the Hipparcos 
Catalogue (while the kinematically improved parallaxes 
from Tycho-2 place the stars at slightly greater distances). 
We have no obvious explanation for this shift, which has 
analogues for other clusters when comparing Hipparcos 
and Tycho-2 solutions. Possibly it reflects the influence of 
subtle systematic effects in the proper-motion data, which 
border on the measurement precision. If this is the case, 
greater confidence should be put on the solution based on 
the Hipparcos data, as the Tycho-2 system of proper mo- 
tions was effectively calibrated onto the Hipparcos system. 
Future space astrometry missions should be able to clarify 
these matters. 

It follows from Eq. ([!]) that any possible bias in the 
estimated space velocity Vq has only a small influence on 
the relative astrometric radial velocities in a given cluster, 
if its angular extent is not too large. Thus the astrophys- 
ical differences would still show up as systematic trends 
when the astrometric radial velocities are compared with 
spectroscopic values. In such a comparison, the bias in 
space velocity would mainly introduce a displacement of 
the zero-point, as mentioned above. 
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Fig. 4. The Hyades: Differences between spectroscopic 
radial-velocity values from the literature, and current as- 
trometric determinations. Systematic differences depend 
on spectral type, and on the projected stellar rotational 
velocity Vsmi. (These dependences are correlated since 
rapid rotation dominates for early-type stars.) An in- 
creased blueshift of spectral lines in stars somewhat hotter 
than the Sun (B — V ~ 0.3-0.5) is theoretically expected 
due to their more vigorous surface convection, causing 
greater convective blueshifts. Gravitational redshifts of 
white-dwarf spectra place them far off main-sequence 
stars. The error bars show the combined spectroscopic and 
astrometric errors; stars with errors > 3 km s _1 are omit- 
ted (except for white dwarfs), as are stars that were not 
retained by kinematic solutions from both Hipparcos and 
Tycho-2 data. 



4.2.1. Hyades: Comparison with spectroscopic data 

The astrophysical potential of astrometric radial veloc- 
ities begins to appear when the accuracy is sufficiently 
high to detect differences relative to spectroscopic val- 
ues (e.g., Dravins et al. 



1999a: Dravins 2001, and refer- 



ences therein). Such differences are expected due to stel- 
lar surface convection ('granulation'): most photons from 
a stellar surface are emitted by hot and rising (thus lo- 
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cally blueshifted) convcctivc elements, which contribute a 
greater number of photons than the cool, dark and sink- 
ing areas. The resulting statistical bias causes a convec- 
tive blueshift, theoretically expected to range from some 
0.2 km s _1 in red dwarfs, ~ 0.4 km s^ 1 in the Sun, to 
1.0 km s" 1 in F-type stars with their more vigorous sur- 
face convection, the precise amount varying among dif- 
ferent spectral lines with dissimilar conditions of forma- 
tion. Gravitational redshifts are expected to vary greatly 
between giants (< 0.1 km s _1 ), main-sequence stars 
(0.5-1 km s _1 ), and white dwarfs (perhaps 30 km s _1 ), 
with almost identical redshifts throughout the spectra. 
Additional effects enter for pulsating stars, stars with ex- 
panding atmospheres, and such with other spectral com- 
plexities. 

Our current accuracies permit such studies to be made 
for the Hyades, and perhaps marginally for a few other 
clusters. In Fig. 4 we show the difference between astro- 
metric radial velocities and spectroscopic measurements 
from the literature. The latter v alues are taken from the 
compilation by Perryman et al. (1998), which mostly are 
measurements by Griffin et al. (1985). The values used 



here are their original measurements, i.e. not applying any 
zero-point or other spectral-type dependent 'corrections' 
(such as were later applied in a convergent-point a nalysi s 
for the Hyades from the same data by Gunn et al. 1988Q . 
The plot also includes white dwarfs, whose astrometric ve- 
locities are from Eq. ([!]) , and where the spectroscopic data 
are 'weighted values' (Ha weighted with twice the weight 
of H/3) from Reid ( p9€j ). 

The errors for Fig. 4 were calculated as the quadratic 
sum of the spectroscopic and astrometric uncertainties, 
where Eq. (|3|) was used for the latter. Most of the errors 
are due to the spectroscopic measurements, and it can be 
noted how the scatter is greater for the [suspected] binary 
stars. 

Over much of the main sequence, convective blue-shifts 
and gravitational redshifts partly cancel one another: an 
increased convective blueshift in hotter stars is partly bal- 
anced by an increased gravitational redshift in these more 
massive stars. Nevertheless, it is theoretically expected 
that the strong increase in the vigour of surface convec- 
tion for middle F-stars (B — V ~ 0.4) should blueshift 
their spectra by ~ 1 km s _1 relative to those of later-type 
G or K-type stars (B — V ~ 1.0). For yet earlier-type 
stars, there do not yet exist any detailed theoretical mod- 
els in the literature from which the convective shift can be 
reliably predicted. 

We believe the expected effects are visible in Fig. 4. 
There is clearly a gradient in the relevant spectral range 
(B — V between 0.4 and 0.7), with roughly the theoreti- 
cally expected sign and magnitude of the effect. The trend 
seems to continue towards even earlier types. 

This is not the first time spectral-type dependent 
radial velocities are seen: a trend of increased spectral 
blueshift in earlier-type stars was already suggested from 
residuals in the convergent-point solution by Gunn et al. 



ants and dwarfs, suggesting differences in gravitational 
redshift, was noted from velocity histograms of giants 
and dwarfs, respectively, in the open cluster NGC3680 by 
Nordstrom et al. fll997| ). 

Both of these works raise an important point relat- 
ing to the sample selection. Spectroscopic velocities are 
usually important for the determination of membership 
probabilities, which are therefore in principle affected by 
systematics of the kind shown in Fig. 4. Spectral shifts 
should therefore be taken into account, lest they influ- 
ence the membership determination and hence the final 
result, including the spectral shifts themselves. Our ini- 
t ial H yades sample is based on that of Perryman et al. 
( 1998 ) , who used their compilation of spectroscopic radial 
velocities to compute membership probabilities. Given the 
relatively large spectroscopic uncertainties for the early- 
type stars, this effect probably did not affect the present 
Hyades sample. However, as long as the mean spectral 
shifts remain unknown, e.g. as a function of spectral type 
along the main sequence, it would be necessary to down- 
weight the more precise spectroscopic velocities in order 
to avoid possible selection effects related to the spectral 
shifts. 

The errors in the spectroscopic velocities in several of 
the hottest (and often rapidly rotating) stars are large, 
making conclusions in that part of the diagram difficult. 
For such stars with often complex spectra and perhaps 
expanding atmospheres, the concept of spectroscopic ra- 
dial velocity must be precisely defined, if studies on the 
sub-km s _1 are to be feasible (cf. Andersen & Nordstrom 
Griffin et al. 



1983 



2000) 



(1988). A difference between the wavelength scale of gi- 



4.2.2. Hyades: The Hertzsprung-Russell diagram 

Already the trigonometric parallaxes from Hipparcos yield 
quite accurate absolute magnitudes, enabling a precise 
Hertzsprung-Russell diagram to be constructed. Our kine- 
matically improved parallaxes permit this to be carried 
further, also verifying the working of our mathematical 
methods. 

Figure 5 shows the gradual improvements in the defini- 
tion of the Hyades main sequence with successively better 
data. The top frame shows the apparent magnitudes (i.e., 
effectively placing all stars at the same mean distance) ; the 
second frame shows the improvement from Hipparcos hav- 
ing been able to resolve the depth of the cluster; the third 
frame uses our kinematically improved parallaxes which, 
especially for the fainter stars, significantly improve the 
definition of the main sequence. This is further marginally 
improved by the use of Tycho-2 data in the bottom frame. 
The plot only shows those stars that were retained in both 
the Hipparcos and Tycho-2 solutions, and excludes [sus- 
pected] binary stars as defined in Sect. ||. 

Besides permitting searches for fine structure in the 
HR diagram, this also confirms the validity of the kine- 
matic solution for the radial- velocity determinations: since 
no photometric information was used in the solution, such 
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Fig. 5. The Hyades: Improved definition of the 
Hcrtzsprung-Russell diagram from kinematically im- 
proved parallaxes. From top: (a) Apparent magnitudes as 
measured by Hipparcos; (b) Absolute magnitudes from 
Hipparcos parallaxes tighten the main sequence since 
the cluster depth is resolved; (c) Absolute magnitudes 
from kinematically improved parallaxes computed from 
Hipparcos data greatly improve the definition of the 
main sequence (especially for the fainter stars), further 
marginally improved by the use of Tycho-2 data (d). This 
permits searches for fine structure in the HR-diagram, 
and also confirms the validity of the kinematic solution 
for the radial-velocity determinations. Only single stars, 
retained by the kinematic solutions from both Hipparcos 
and Tycho-2 data, are plotted. 



an improvement could hardly be possible unless the un- 
derlying physical model is sound. With the kinematically 
improved parallaxes the error in Mh p is typically only 
~ 0.03 mag, much smaller than the symbol size in Fig. ^[ 

In addition to the two giants retained in the solutions, 
a few stars lie off the main sequence. Below it is HIP 10672 
at B — V = 0.567, a single star quite far (some 30 pc) from 
the cluster centre; HIP 17962 at B -V = 0.782, an eclips- 
ing binary containing a hot white dwarf (Nelson & Young 



1970) which causes a displacement towards the blue; and 
HIP 19862 at B - V = 0.924, with an uncertain colour 
index in the Hipparcos Catalogue [gb-v = 0.301) - the 
value B — V = 1.281 ± 0. 009 g iven in the Hipparcos Input 



Catalogue (Turon et al. 1992| ) would place it exactly on 
the main sequence. All these stars have an uncertainty in 
Map of 0.05 mag or less, meaning that they are not mis- 
placed vertically. Simulations of a Hyades-type cluster by 
Portegies Zwart et al. ( 2001 ) showed that a few stars end 
up below the main sequence as a result of binary inter- 
action leading to blueward displacements. From Fig. 5 it 
is difficult to tell where the turnoff point really is: some 
stars to the far left may be blue stragglers. 

The rest of the cluster stars lie practically on a sin- 
gle curve, which can be considered a confirmation of the 
parallax improvement. It is not clear whether the remain- 
ing spread of the main sequence in Mu P is real or can be 
accounted for by uncertainties in B — V , although these 
are small. Effects such as differential reddening within the 
cluster seem unlikely: Taylor (1980) found only a very 
small colour excess E(B - V) = 0.003 ± 0.002 mag for 
the Hyades. 

Improved absolut e ma gnitudes were also determined 
by de Bruijne et al. ( 2001 ) based on our o riginal method 
(Dravins et al. 1997 ). Lebreton et al. ( 2001 ) compared our 
kinematically improved parallaxes to those by de Bruijne 
et al., finding excellent agreement in all values, except for 
one star (HIP 28356). We note that this particular star is 
the one located the furthest from the cluster centre, and is 
also one where our cleaning procedure removed it from the 
Tycho-2 solution (although it was retained in Hipparcos 
data) . It may be a long-period binary whose photocentric 
motion causes a deviation in the modulus of the measured 
proper motion, if not in its direction. 

For further discussions of the post-Hipparcos HR dia- 



gram for the Hyades, see Perryman et al. (1998), Madsen 
( |1999| ), Castell ani et al. ( |200lD , de Bruijne et al. ( |200lD , 
and Lebreton ( pOOl ), where the latter four have used the 
improved parallaxes. 



up 



4.3. The Coma Berenices cluster 

The Coma Berenices sample is made 
Hipparcos stars in Odenkirchen et al. 
pie includes four stars that, while slightly beyond their 
selected limit for membership, nonetheless were consid- 
ered to 'very probably also belong to the cluster'. Since 
the small number of stars made the solution unstable al- 



of the 40 
(|1998|). This sam- 
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ready after rejecting two of them, the results in Table [j] 
are given for the full sample (gu m — oo). Although the 
typical errors in the astrometric radial velocities are only 
1.2 km s -1 , the precision of published spectroscopic values 
is generally insufficient for meaningful comparisons. 

The kinematically improved parallaxes produce only a 
slight improvement in the HR diagram at the red end of 
the main sequence (Fig. ^|). 

4.4. The Pleiades cluster 



The sample contains 60 stars from van Leeuwen ( 1999 , and 
private comm.). The stars are too few and/or the cluster 
too distant for the basic cluster model to give a direct 
estimate for the velocity dispersion. It has instead been 
estimated by the procedure described in Appendix A. 4 of 
Paper II. No improvements to the parallaxes result from 
the kinematic solution, since our method is unable to re- 
solve the depth of this cluster; it therefore in essence as- 
cribes the same distance to every star. The small angular 
extent means that also the astrometric radial velocity is 
practically the same for all stars (Fig. ||[ bottom) . 

As seen in Fig. ^, the Pleiades main sequence occupies 
a position at the lower edge of the distribution for the 
different clusters. The Pleiades cluster is at the focus of an 
ongoing debate concerning possible localized systematic 
err ors in the H ipparcos paralla xes (se e e.g. Pinsonneault et 



al. 1998, 200C; Robichon et al. 1999: Narayanan & Gould 



199S; van Leeuwen 199S, 2000; Paper II; Stcllo & Nissen 



2001 ) . If such errors were present in our input data, they 



would not be detected by the present maximum-likelihood 
method, but would affect also the kinematically improved 
parallaxes. Consequently, the present results provide no 
direct new information towards the resolution of this issue. 



4.5. The Praesepe cluster 

The inves tigate d sample was based on 24 stars from van 
Leeuwen (199E, and private comm.). As for the preceding 
two clusters, the solution places the stars at practically 
the same distance. The resulting mean astrometric radial 
velocity has an error of ~ 15 km s _1 . While demonstrating 
the applicability of the method, this present accuracy is 
insufficient for detailed stellar studies. 



4.6. A composite HR diagram 

Kinematically improved parallaxes from the different clus- 
ters enable a very detailed comparison between the main 
sequences of different clusters. Such a Hertzsprung-Russell 
diagram for five nearby clusters is in Fig. ^. We again 
stress that, while our kinematic solution reduces the ran- 
dom noise, it does not address any possible systematic 
effects and therefore cannot decide whether, e.g. the sys- 
tematic shifts in luminosity between different clusters are 
caused by astrophysical or by instrumental effects. 
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Fig. 6. Hertzsprung-Russell diagram for single stars in the 
better-defined open clusters, obtained using kinematically 
improved parallaxes from Hipparcos data. Hp magnitudes 
are given since these are more precise than ground-based 
V ones, and since Hp values are available for all stars 
in this sample (values for M# p fall quite close to those 
of My). The random errors in these kinematically im- 
proved parallaxes are lower by factors of typically 2 or 3 
compared with the original Hipparcos values, and the ab- 
solute magnitudes are correspondingly more precise, be- 
ginning to reveal fine structures in the HR diagram. For 
the Hyades, de Bruijne et al. ( 2000| , 2001 ) suggested the 
existence of two underpopulated main-sequence segments 
around B — V ~ 0.38 and 0.5, identified as 'Bohm-Vitense 
gaps', theoretically predicted due to changing efficiencies 
of stellar convection at the corresponding temperatures. 
However, these gaps are not seen in other clusters, and 
their 'existence' is consistent with small-number statistics 
causing random clustering. This probably also applies to 
the apparent 'Gap 3' at B — V ~ 0.8. 



However, the very low noise level permits to search for 
morphological fine structures in the HR diagram. From 
post-Hippar cos data for the Hyades, de Bruijne et al. 
( |200uj 2001 ) suggested the existence of two underpopu- 
lated main-sequence segments around B — V approx 0.38 
and 0.5, identified as 'Bohm-Vitense gaps', theoretically 
predicted due to changing efficiencies of stellar convection 
at temperatures corresponding to those particular colours. 
However, these gaps are not seen in other clusters, and 
their 'existence' is consistent with small-number statistics 
causing random clustering. This probably also applies to 
the apparent 'Gap 3' at B — V ~ 0.8. (Of course, the 
detectability of such gaps depends also on the precision 
in the other axis, i.e. the colour index). Although some 
authors have suggested a possible presence of such 'gaps' 
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(e.g., Rachford & Canterna 2000), extensive analyses of 
field stars, using Hipparcos parallaxes, failed to show any 
(Newberg & Yanny |1998fl . 

A certain fine structure (wiggles, etc.) is theoretically 



expected in the HR- diagram (e.g., Siess et al. 1997); pos- 
sible hints of that are becoming visible for the Hyades. 

5. Application to OB associations 

The procedures of determining astrometric radial veloci- 
ties were applied also to a number of nearby associations 
of young stars. The situation is here somewhat different 
from that of the previously discussed older clusters be- 
cause (at least some of) these younger associations may be 
undergoing significant expansion, or have otherwise com- 
plex patterns of stellar motion on levels comparable to 
our desired accuracies. We recall that the present moving- 
cluster method is based upon measuring the rate of angu- 
lar expansion or contraction: it cannot therefore segregate 
whether a change in angular scale occurs because the clus- 
ter is approaching or expanding. While - on the accuracy 
levels aimed at - this should not be a problem for the 
older clusters, the likely expansion of young associations 
may introduce significant biases in the solution. Another 
complication is that, since some of the associations cover 
large areas of sky, there is an increased risk for contamina- 
tion of the samples by field stars. Further, spectroscopic 
radial velocities often cannot be used to decide member- 
ship, both because they do not exist in significant num- 
bers, and because their actual measurement is difficult for 
the often rapidly rotating O and B-type stars that make 
up much of these associations. For such reasons, the asso- 
ciations are here being treated separately. 

The predictions in Paper I indicate that the accura- 
cies of Hipparcos should enable astrophysically interest- 
ing results to be obtained for perhaps half a dozen of the 
nearer associations. Among these, Lower Centaurus Crux, 
Upper Centaurus Lupus and Upper Scorpius form part of 
the larger Scorpius OB2 complex, while the a Persei and 
'HIP 98321' associations are indepe nden t entities. 

Except for 'HIP 98321' (Sect. |1§ the selection of 
members in the different associations is based on data 



from de Zeeuw et al. ( 1999 ). In their sample selection they 
combine one method using Hipparcos positions and proper 
motions (de Bruijne 1999a ), and another using Hipparcos 
positions, proper motions and parallaxes (Hoogerwerf & 
Aguilar 1999 ). Although this could cause some contamina- 
tion by outliers, simulations showed that only 20% of the 
stars in the first method are expected to be field stars, and 
only 4% in the second. Although, in principle, our proce- 
dure for rejecting outliers does reduce this contamination, 
actually only few stars were rejected. 

De Bruijne ( 1999b] ) used an implementation of our 
original method (Dravins et al. 1997) to obtain kinemati- 
cally improved parallaxes for the three OB associations in 



some internal structure. We refer to de Bruijne's work 
concerning the three-dimensional structure of the com- 
plex, although his distance estimates are slightly different 
from ours (mainly because his selection criterion, <7i; m = 9, 
differs from our gn m — 15). Based on the kinematically im- 
proved parallaxes, we presented the Hertzsprung-Russell 
diagrams of Upper Centaurus Lupus and Lower Centaurus 
Crux in Madsen et al. (200C). For additional discussion 
of the HR diagrams of the complex we again refer to de 
Bruijne ( |l999b| ). 

The solutions for the associations as a whole were given 
in Table [l] (and its electronic version) , while the results for 
the individual stars are given in the electronic version of 
Table 0. 



5.1. The Lower Centaurus Crux association 

The cleaned sample consists of 179 stars with an esti- 
mated internal dispersion of 1.1 km s . Combined with 
the rather small uncertainty of the cluster velocity, the re- 
sulting standard error for the astrometric radial velocities 
is 1.2 to 1.3 km s _1 . In the corresponding HR diagram, the 
main sequence becomes somewhat better defined, most 
noticeably in the A-star regime (Madsen et al. |2000D , al- 
though there still remains a significant spread. 



5.2. The Upper Centaurus Lupus association 

For Upper Centaurus Lupus, the rejection procedure pro- 
duces a clean sample with 218 stars with an estimated 
internal dispersion of 1.2 km s _1 . The resulting standard 
error of the astrometric radial velocities is 1.3 km s . 

The HR diagram clearly shows an improvement across 
the whole spectral range of the main sequence (Madsen 
et al. 200C). Probably, the remaining spread is caused 



by non-detected binaries, some non-members, and pre- 
main sequence objects moving onto the main sequence. 
Differential reddening across the association and perhaps 
also in depth could also cause a spread of the main se- 
quence, although Upper Centaurus Lupus is not believed 
to be as much affected as the other two associ ations in the 
Scorpius OB2 complex (de Zeeuw et al. 1999] ). 



the Scorpius OB2 complex (cf. Sect. 3.6). While the depth 
of the associations is not fully resolved by the Hipparcos 
parallaxes, the kinematically improved parallaxes reveal 



5.3. The Upper Scorpius association 

The maximum-likelihood solution for the Upper Scorpius 
association became unstable after rejection of eight stars, 
at which point the criterion g max < 15 was still not met. 
We therefore choose to give results for the solution using 
all 120 stars in the original sample. The internal dispersion 
is in line with that of the previous two associations, but 
the larger uncertainty in the cluster velocity gives a higher 
standard error of about 1.9 km s _1 for the astrometric 
radial velocities. 

It is difficult to judge whether the main sequence is 
actually better delineated by the kinematically improved 
parallaxes. Upper Scorpius appears to be close to the limit 
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Table 2. Data for individual stars for the better-defined clusters. Estimated radial velocities and their standard errors 
are derived from Eqs. (||) and (||), using the adopted solutions in Table |l| and the corresponding covariances. Columns: 
HIP = Hipparcos Catalogue number, v r ^i p = astrometric radial velocity [km s _1 ] obtained from the kinematic solution 
using data from the Hipparcos main catalogue; 7?Hip = estimated parallax [mas] from the kinematic solution using 
data from the Hipparcos main catalogue; e(7?Hip) = standard error [mas] of this estimated parallax. Part 1: data for 
Ursa Major and Hyades. The complete table for all clusters and associations, including results and errors obtained 
from both Hipparcos and Tycho-2 data, is available in electronic form. 



HIP iv.Hip *Hip e(*Hip) HIP « r ,Hip *Hip e(*Hip) HIP « r ,Hip *Hip <s(*Hip) HIP « r> Hip *Hip <s(*Hip) 



Ursa Major 



2213 


3 


10 


6 


6 





6 


38228 


-17 


73 


45 


8 





9 


61100 


—12 


17 


41 


1 


2 


9 


77233 


2 


01 


21 


3 


(1.8 


8486 


3 


75 


42 





3 


8 


42438 


-15 


45 


70 








7 


61481 


-11 


96 


38 


1 


1 


1 


80337 


12 


39 


77 


6 


(1.9 


8497 


3 


74 


42 


4 





9 


43352 


-6 


85 


14 





1 


1 


61621 


-1 


0(1 


20 








7 


80686 


12 


20 


82 





(1.0 


10403 


-7 


28 


24 


4 





9 


46298 


-12 


80 


19 


1 





9 


61946 


-11 


83 


43 








8 


80902 


-8 


57 


30 


7 


(1.0 


10552 


-0 


98 


28 


9 


2 


7 


48341 


-12 


85 


10 


1 





8 


62512 


-11 


05 


40 


8 


3 





82780 


4 


(10 


7 


2 


(1.0 


17874 


1 


39 


15 


8 





6 


48356 


-10 


57 


11 








8 


62956 


-11 


42 


40 


3 





6 


83988 


-2 


42 


40 


5 


1.8 


18512 


—5 


98 


63 


3 


2 





49593 


-16 


04 


35 


7 





8 


63008 


-8 


70 


28 


3 


1 


5 


83996 


-2 


42 


40 


1 


2.8 


19655 


-11 


96 


38 


5 


1 


3 


49929 


-15 


06 


15 


1 





8 


63503 


-11 


21 


40 








6 


87079 


-6 


33 


25 


1 


(1.0 


10855 


_g 


28 


47 


g 








—15 


88 


12 






g 


04405 


g 


01 


]4 






2 




10 


(12 


57 




(1 8 


19859 


-8 


28 


47 


1 


1 


1 


51814 


-15 


21 


37 


8 





6 


64532 


-10 


78 


38 








7 


91159 


9 


81 


28 


9 


1.5 


21818 


-11 


0!) 


74 


6 


1 


2 


53910 


—14 


08 


41 


1 








65327 


-10 


50 


39 


8 


1 


4 


94083 


-6 


00 


36 





0.5 


23875 


-9 


03 


30 


7 





8 


53985 


-14 


19 


85 


7 


1 


4 


65378 


-10 


34 


41 


7 








96258 


1 


18 


39 


2 


0.5 


25110 


-11 


98 


47 


7 





5 


55454 


-6 


16 


75 


8 


1 


7 


65477 


-10 


30 


40 


2 








101027 


17 


21 


33 


(J 


(1.9 


27072 


-6 


83 


111 


5 





6 


56154 


-7 


38 


15 


7 





8 


66459 


-8 


0!) 


01 


7 


1 


2 


103738 


17 


75 


14 


5 


0.8 


27913 


-14 


70 


115 


4 


1 


1 


57283 


-5 


56 


9 


4 





7 


69989 


-2 


20 


38 


2 





8 


106481 


2 


56 


20 


2 


0.5 


28954 


-14 


70 


04 


7 





9 


57548 


-9 


22 


209 


9 


2 


2 


71876 


-8 


19 


23 


5 








110091 


12 


91 
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1 


(1.9 
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-5 


20 


13 
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6 
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-13 


38 


30 








7 
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61 
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7 


1 


7 
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37 
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2.(1 


30630 
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68 


3 


1 
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-12 


73 


35 


3 


1 


2 
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75312 




58 


5 3 
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40 
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22 


8 
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32 
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3 
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38 


03 


18 
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2 
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71 


21 
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3 
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41 
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20 


4 


0.4 
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25 
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(i 


20130 


38 


34 


21 


9 


(1 


3 


20826 


39 


99 


22 
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4 


21683 


40 


77 


18 


3 


0.3 


15563 


31 


72 


32 


1 





4 


20146 


38 


00 


21 
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4 


20827 


39 


82 


20 


5 





4 


21723 


41 


08 


22 


9 


0.0 


15720 


30 


(i() 


31 








5 


20205 


38 


91 
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1 





3 


20842 


38 
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20 
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39 


75 


16 


6 
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03 
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07 
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20349 


38 


43 


20 


2 





3 


20899 


39 


05 


21 








3 


22253 


40 


40 


18 


3 


0.4 


18170 


35 


77 


23 
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3 


20350 


38 


80 


21 


4 





3 


20901 


40 


02 


21 


3 





3 


22265 


41 


21 


20 





0.4 


18322 


30 


32 


21 


5 





4 


20357 


39 
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20 
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3 


20916 


39 


79 


18 


7 





5 


22271 


39 


70 


27 


3 


0.5 


18327 


36 


03 


24 








20400 


30 


27 


22 


3 




3 


20035 


30 


07 


2 1 


g 




3 


99^0 


40 


80 


9(1 


g 


(J 4 


18658 


36 


94 


23 


5 
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20419 


39 


40 


22 


1 


(1 
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20948 
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05 


21 


8 


(1 
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22380 


41 


27 


20 


9 


0.4 
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36 


58 
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20440 
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27 
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5 


20949 


38 


12 
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3 


22394 


40 


18 


20 


4 


0.4 


18946 


36 


70 


21 








4 


20455 


39 


04 


21 


1 





3 


20951 


39 


65 


22 


3 





4 


22422 


41 


05 


20 


8 


0.4 


19082 


30 


96 


20 


8 





5 


20480 


38 


50 


10 


8 





3 


20978 


39 


83 


22 








4 


22505 


41 


81 


21 


9 


0.4 


19098 


37 


10 


22 


1 





4 


20482 


38 


82 


10 


2 





3 


20995 


39 


95 


22 


3 





4 


22524 


41 


72 


20 


3 


0.4 


19148 


37 


44 


20 


9 





3 


20484 


39 


17 


20 


8 





3 


21008 


39 


40 


10 


1 





3 


22550 


42 


15 


21 


3 


0.4 


19207 


37 


56 


21 


6 





4 


20485 


39 


27 


24 


9 


(1 


5 


21029 


39 


93 


21 


8 





3 


22565 


41 


43 


19 


3 


0.3 


19261 


37 
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22 
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3 


20491 


38 


04 


18 


7 





3 


21036 


40 


15 


22 


4 





3 


22566 


41 


85 


10 


4 


0.3 


19263 


37 


53 


21 


7 





4 


20492 


39 


38 


21 








4 


21039 


39 


99 


21 


7 





4 


22654 


41 


48 


19 


2 


0.5 


19316 


37 


91 


20 
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5 


20527 


39 


40 


22 


6 
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21066 
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34 


21 
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41 
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35 


57 
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39 
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15 
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12 
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5 


21123 


39 


87 


22 
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4 


23312 


42 


94 


19 


3 


0.5 
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38 


28 


20 


5 





4 


20567 


39 


24 


10 


8 





4 


21137 
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09 


22 
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41 


84 
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0.3 
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8 
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21152 
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00 
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3 


21256 


39 


50 


23 
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4 


23983 


43 


56 


19 





0.4 
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37 


31 


22 


4 





3 


20641 


38 


02 


22 


5 





3 


21261 


39 


89 


21 


6 





5 


24019 


40 


81 


17 


9 


0.3 


19796 


38 


65 


22 








3 


20648 


39 


25 


21 


8 





3 


21267 


40 


49 


21 


8 





4 


24116 


42 


53 


12 


4 


0.3 


19808 


38 


58 


21 


9 
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20661 


39 


48 


21 


3 


(1 


3 


21273 


40 


37 


22 


4 





5 


24923 
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19834 


38 


53 


21 


3 





8 
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39 


27 


22 


9 
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21317 


40 


38 
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26382 


44 


49 
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7 


0.5 
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38 


40 


21 


8 
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39 


17 


22 


1 





4 


21459 


39 


43 


23 


6 





3 


28356 


45 


83 


13 


4 


0.0 



of our method, due to its larger distance, smaller angular 
size and a smaller number of member stars, compared with 
Lower Centaurus Crux and Upper Centaurus Lupus. 



5.4. The Scorpius OB2 complex 

Lower Centaurus Crux, Upper Centaurus Lupus and 
Upper Scorpius are all part of a larger OB complex, 
known as Scorpius OB2, with similar space velocity vec- 
tors (Blaauw 1964). Therefore, an attempt was also made 
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Table 2. (Continued) Data for individual stars for the better-defined clusters. Part 2: data for Coma Berenices, 
Pleiades and Praesepe. The complete table for all clusters and associations, including results and errors obtained from 
both Hipparcos and Tycho-2 data, is available in electronic form. 



HIP 


Hip 


■""Hip 


e(*Hip) 


HIP 


"r.Hip 


^Hip 


f(*Hip) 


HIP 


"r.Hip 


"^Hip 


s(*Hip) 


HIP 


"r.Hip 


fillip 


e(* H 


Coma Berenices 




























59364 


-1.35 


11.0 


0.6 


60123 


-1.59 


11.5 


().(i 


60525 


-1.48 


10.7 


O.fi 


61205 


-2.25 


13.2 


1.0 


59399 


-1.27 


10.5 


0.8 


60206 


-1.43 


11.3 


0.7 


60582 


-1.49 


10.3 


0.7 


61295 


-1.56 


10.4 


0.5 


59527 


-1.39 


11.4 


0.6 


60266 


-1.39 


11.0 


0.0 


60611 


-1.34 


11.4 


0.7 


61402 


-1.79 


11.0 


0.7 


59833 


-1.34 


10.4 


0.6 


60293 


-1.50 


10.9 


0.8 


60649 


-1.59 


12.0 


0.6 


62384 


-1.70 


11.9 


(1.7 


59957 


-1.36 


11.3 


0.6 


60304 


-1.55 


11.1 


0.7 


60697 


-1.63 


13.1 


0.5 


62763 


-2.10 


10.9 


0.0 


60014 


-1.85 


12.3 


0.6 


60347 


-1.23 


11.2 


0.6 


60746 


-1.61 


11.6 


0.5 


62805 


-1.98 


12.5 


0.8 


60025 


-1.21 


12.2 


0.7 


60351 


-1.46 


11.4 


0.5 


60797 


-1.56 


10.9 


0.6 


63493 


-2.02 


11.7 


0.7 


60063 


-1.33 


11.5 


0.7 


60406 


-1.48 


10.6 


(J. 7 


61071 


-1.53 


11.2 


0.5 


64235 


-2.78 


12.8 


1.(1 


60066 


-1.41 


12.3 


0.5 


60458 


-1.56 


11.6 


0.6 


61074 


-1.74 


11.9 


0.6 


65466 


-2.48 


11.7 


0.6 


60087 


-1.21 


11.6 


0.5 


60490 


-1.49 


11.8 


0.5 


61147 


-1.86 


11.6 


0.6 


65508 


-2.60 


11.2 


0.8 


Pleiades 






























16217 


9.06 


8.6 


0.2 


17289 


11.20 


7.9 


0.3 


17573 


10.70 


8.5 


0.2 


17847 


10.98 


8.3 


(1.2 


16407 


9.63 


8.7 


0.2 


17317 


11.40 


8.4 


0.3 


17579 


10.62 


8.4 


(1.2 


17851 


10.95 


8.7 


(1.2 


16635 


9.93 


8.2 


0.4 


17325 


12.39 


8.7 


0.2 


17583 


10.25 


8.5 


0.2 


17862 


10.85 


8.3 


0.2 


16639 


10.17 


8.2 


0.3 


17401 


10.92 


8.5 


0.2 


17588 


10.64 


8.3 


0.2 


17892 


11.67 


8.4 


0.2 


16753 


10.90 


8.7 


0.3 


17481 


12.23 


8.3 


0.2 


17607 


12.25 


7.9 


0.4 


17900 


11.16 


8.2 


0.2 


16979 


11.18 


8.1 


0.3 


17489 


10.69 


8.3 


0.2 


17608 


10.91 


8,1 


(1.2 


17923 


11.11 


8.2 


0.3 


17000 


11.12 


8.4 


0.2 


17497 


11.14 


8.4 


(J. 2 


17625 


10.07 


8.4 


0.2 


17999 


11.10 


8.4 


(1.2 


17020 


10.30 


8.1 


0.3 


17499 


10.77 


8.5 


0.2 


17664 


10.68 


8.6 


0.2 


18050 


10.87 


8.5 


0.2 


17034 


10.25 


8.6 


0.2 


17511 


11.69 


8.1 


(J. 3 


17692 


11.01 


8.3 


0.2 


18154 


10.84 


8.3 


(1.3 


17044 


10.36 


8.4 


0.3 


17525 


9.81 


8.9 


0.5 


17694 


11.40 


8.5 


0.3 


18263 


11.07 


9.1 


0.4 


17091 


10.90 


8.8 


0.3 


17527 


10.46 


8.6 


0.2 


17702 


10.89 


8,1 


0.2 


18266 


12.38 


8.3 


0.4 


17125 


8.93 


8.1 


3 


17531 


10 63 


7 9 


(j 2 


17704 


10 81 


8 3 


2 


18431 


11 69 


8 (i 


(1 2 


17168 


8.49 


7.5 


0.4 


17547 


8.75 


8.9 


0.3 


17729 


10.33 


8.6 


(1.2 


18544 


12.83 


9.5 


0.3 


17225 


10.89 


8.6 


0.3 


17552 


12.05 


9.0 


0.2 


17776 


11.23 


8.7 


(1.2 


18559 


11.35 


8.2 


0.2 


17245 


9.95 


8.0 


0.3 


17572 


11.24 


8.5 


0.2 


17791 


10.82 


8.0 


(1.2 


18955 


12.06 


8.5 


0.3 


Praesepe 






























41788 


37.62 


5.9 


0.3 


42327 


36.63 


5.7 


0.3 


42556 


36.25 


5.5 


0.2 


42952 


35.76 


5.3 


0.3 


42133 


36.86 


5.7 


0.3 


42485 


36.30 


5.7 


0.3 


42578 


36.19 


5.3 


(1.2 


42966 


35.47 


(i.O 


0.3 


42164 


36.78 


5.7 


0.3 


42516 


36.21 


5.4 


0.2 


42600 
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to combine the three associations in a single solution as- 
suming a common velocity vector. The resulting vector 
and internal dispersion are in Table |, but we give no re- 
sults for individual stars. 

The HR diagram is visibly improved, indicating that 
Sco OB2 could meaningfully be regarded as one single 
structure. However, a combination of the HR diagrams 
from the three separate solutions is even slightly better 
defined, suggesting that Sco OB2 is, after all, better con- 
sidered as three separate structures. 

When treating Sco OB2 as one complex, the estimated 
internal velocity dispersion is only slightly larger than for 
the separate solutions, and the formal uncertainty of the 
space velocity vector is remarkably small. Nevertheless, 
when comparing the resulting astrometric radial veloci- 
ties with those from the previous solutions we find notice- 
able differences. For Lower Centaurus Crux (LCC) we find 
(u„(LCC) - v„(Sco OB2)) ~ -2 km s" 1 , while for Upper 
Centaurus Lupus and Upper Scorpius the corresponding 
mean differences are +4 km s _1 and +10 km s _1 , respec- 
tively. Such a progression of systematic differences could 
be expected if Sco OB2 is not really one uniform complex, 
or if there is some internal velocity field. At any rate, the 
comparison shows that one has to be careful when inter- 
preting the results for young associations: although we get 
a stable solution with small residuals when considering the 



whole complex, the resulting velocities are not trustwor- 
thy. 

Thus both the HR diagrams and the radial-velocity 
solutions indicate that the Sco OB2 complex has some in- 
ternal kinematic structure that ultimately will need to be 
modelled, although it is only marginally discernible in the 



present data. In Sect. 5.7 we discuss the possible expan- 
sion of the associations. 



5.5. The a Persei association (Per OB3) 

This a Per association is sometimes denoted an open clus- 
ter. From our sample we obtain a mean astrometric ra- 
dial velocity of 4.5 ± 2.2 km s^ 1 . A rather modest inter- 
nal velocity dispersion a v ~ 0.7 km s _1 was found using 
the procedure of Appendix A. 4 in Paper II. The value 
is smaller than for the other OB associations, and indi- 
cates that it may be reasonable to look upon the struc- 
ture as a young open cluster instead. The velocity dis- 
persion, together with the uncertainty in the solution for 
the cluster velocity, combine to give a standard error of 
about 2.3 km s -1 in the astrometric radial velocities of the 
individual stars. The parallax improvement is not good 
enough to have a visible impact on the HR diagram. 

Our radial-velocity result is close to the spectroscopic 



values of ~ +2 km s 1 (Prosser 1992), while somewhat 
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larger than the —0.9 km s 1 derived from the convergence- 
point solution by Eggen (1998). 



5.6. The 'HIP 98321' association 

This possible association was recently discovered in the 
Cepheus-Cygnus-Lyra-Vulpecula region by Platais et al. 
( 1998 ), during a search for new star clusters from 
Hipparcos data. They named it after the central star 
HIP 98321, and found 59 probable members. Because of 
the Hipparcos limiting magnitude, only O, B, and A-type 
stars are utilizable. It was a bit surprising to find that 
this association gives a good kinematic solution despite 
its great distance of 307 pc; the reason is probably its 
large mean radius on the sky of ~ 12 degrees (Fig. 1). 

The mean astrometric radial velocity is —19.3 ± 
1.6 km s" 1 for the sample of all 59 stars. Together with 
the estimated internal dispersion of 2.6 km s _1 , the stan- 
dard error of the individual astrometric radial veloci- 
ties is around 3.2 km s _1 . These values are consistent 
with the somewhat uncertain spectroscopic velocities for 
these early-type stars. Published values spread around 
— 15 km s , suggesting a possible expansion of the as- 
sociation compatible with its isochrone age (Table 4 in 
Paper I and next section). 

During the cleaning process, the maximum g value was 
always below gn m = 15; thus no star was removed from the 
original sample. Some contamination by outliers may nev- 
ertheless be expected due to the lack of spectroscopic in- 
formation in the selection of the stars. It would have been 
a nice confirmation of the existence of this new associa- 
tion if the improved parallaxes had given a better-defined 
main sequence, but unfortunately the improvement is not 
sufficient to have any visible effect in the HR diagram. 



5.7. Expanding associations? 

Figure [?] shows the astrometric versus spectroscopic ra- 
dial velocities for stars in the Scorpius-Centaurus group 
of young associations, both for each individual subgroup, 
and for the complex treated as a whole. The spectro- 
scopic values are those compiled in the Hipparcos Input 



Catalogue HIC (Turon et al. 1992). Because we are deal 



ing with young and rapidly rotating early-type stars, the 
spectroscopic errors are quite large; some contamination 
is also expected due to outliers and binaries. 

The astrometric radial velocities in OB associations are 
expected to show a significant bias due to expansion ef- 
fects (Paper I). Assuming the inverse age of an association 
to be the upper limit on the relative expansion rate, the 
resulting maximum bias in the astrometric radial velocity 
can be computed from Eq. (10) in Paper I. This effect is 
directly proportional to the distance to the stars in the 
association and inversely proportional to its age. The ex- 
pansion causes a positive shift in iv(spectr) — tv(astrom): 
the cluster's increasing angular size is wrongly interpreted 
as approaching motion. This [upper limit of the] expan- 
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Fig. 7. Astrometric versus spectroscopic radial velocities 
for stars in the Scorpius-Centaurus group of young associ- 
ations, expected to undergo kinematic expansion. The top 
three frames show separate solutions for each subgroup. 
Assuming a rate of isotropic expansion equal to the inverse 
age of the cluster, a bias in the astrometric radial velocity 
would result, marked by dashed gray lines (the cluster's 
increasing angular size would be interpreted as approach- 
ing motion; Paper I). The assumed ages are 11, 14 and 
5 Myr, respectively (de Geus et al. 1989). The bottom 
frame shows the solution for all 510 stars in the groups, 
treated as one entity (only stars with known spectroscopic 
velocities are plotted). While these data do indicate some 
expansion of this complex of young associations, the ex- 
pansion of its individual parts is significantly slower than 
the naively expected rate. 
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sion bias is plotted in Fig. [7] together with the spectro- 
scopic and astrometric velocities. We have not been able 
to observe any correlation between distance and expan- 
sion with the present data. The expected effect should be 
a few km s -1 , but it probably drowns in the noise from 
spectroscopic measurements that have errors of compara- 
ble magnitude, and from a possible anisotropic expansion. 

The interpretation of Fig. |?] is not obvious. Stars in 
Lower Centaurus Crux show a wide spread in the spectro- 
scopic values, while the mean is roughly consistent with an 
isotropic expansion at about half the rate naively expected 
from the age of the association. The same can be said for 
Upper Centaurus Lupus. For these associations the indi- 
cated 'kinematic age' (equal to the inverse of the current 
expansion rate) is thus around 20-30 Myr, or twice the 
isochrone ages according to de Geus et al. (198£). Upper 
Scorpius on the other hand, which is the youngest of the 
subgroups (5-6 Myr according to de Geus et al.), does not 
seem to expand at all: taken at face value, the data rather 
suggest that it contracts. The combined sample again in- 
dicates some expansion, roughly consistent with a kine- 
matic age of 20 Myr. We note that already Blaauw (1964) 
derived such an expansion age of 20 Myr for the Scorpius- 
Centaurus complex as a whole from the ~ 10 km s _1 dis- 
crepancy between the spectroscopic radial velocities and 
the proper motion data combined with photometric dis- 
tances. 

Of the three subgroups in Sco OB2, the result for 
Upper Scorpius thus stands out as rather puzzling. A de- 
tailed study of this association by Prcibisch & Zinnecker 
( |l999l ) suggested that the star formation process was trig- 
gered by a giant supernova explosion in the neighbouring 
Upper Centaurus Lupus. What effect that may have had 
on the internal kinematics of Upper Scorpius is hard to 
say. There is a priori no reason to expect Upper Scorpius 
to be a bound system without expansion. The star forma- 
tion in Upper Scorpius itself seems to have dispersed the 
rest of the parent molecular cloud. This resul t seem s to 
imply the standard picture (see e.g. Mathieu 1986 for a 
review) : the removal of gas leads to loss of binding mass 
of the system, it becomes unbound and consequently will 
expand. 

From calculations inspired by our method, Makarov 
& Fabricius (2001) estimated an expansion rate of 



0.12 km s pc for the TW Hya association of young 
stars, assuming a uniform expansion. The TW Hya associ- 
ation is dominated by late- type stars and may be an exten- 
sion of Lower Centaurus Crux. The expansion corresponds 
to a bias of the centroid radial velocity of ~ —9 km s _1 
- comparable to the biases we find for Upper Centaurus 
Lupus and Lower Centaurus Crux - and to a dynamical 
age of 8.3 Myr, in agreement with previous age determina- 
tions for TW Hya's T Tauri members (Webb et al. |1999| ). 

The results we find here are promising in the sense 
that it is possible to obtain information about the internal 
kinematics, formation history and age, but at the same 
time they confirm the complexity of the kinematics of the 
associations in the Sco OB2 complex. In the end more 



accurate spectroscopic observations are also required to 
answer these questions. These would in particular allow 
true expansion to be disentangled from the perspective 
effects of the radial motion. 



6. Conclusions 

The radial motions of stars have been studied through 



spectroscopy since the year 1868 (Hearnshaw 1986). 
Recently, the accuracies realized in astrometry have en- 
abled such determinations to be made also through purely 
geometric methods. Once sufficient accuracies are reached, 
this will enable an absolute calibration of the stellar ve- 
locity scale for stationary and variable stars, irrespective 
of any complexities in their spectra. Indeed, for several 
early-type stars (with complex spectra smeared by their 
rapid rotation) the radial velocities already now deter- 
mined through astrometry are more accurate than has 
been possible to reach spectroscopically in the past. 

The differences between these astrometric radial- 
velocity values and wavelength measurements of differ- 
ent spectral features may become a new diagnostic tool 
in probing the dynamic structure of stellar atmospheres. 
Already the present work has made available quite accu- 
rate astrometric radial velocities for stars of many more 
spectral types than those for which hydrodynamic model 
atmospheres have been developed (from which, e.g., con- 
vective and gravitational wavelength shifts in their spectra 
could have been predicted). For such stars, the limitations 
in understanding the differences between astrometric and 
spectroscopic radial velocities may now lie primarily with 
spectroscopy and atmospheric modelling, rather than in 
astrometry. 

In this series of three papers, we started by explor- 
ing different types of fundamental possibilities of astro- 
metrically determining radial velocities, identifying which 
methods could be applicable on existing data already to- 
day. Among the latter, the moving-cluster method was 
found capable of yielding astrophysically interesting, sub- 
km s _1 accuracies, and its mathematical methods were 
developed in Paper II. In the present paper, data from 
Hipparcos were used in applying the method to obtain 
solutions for more than 1 000 stars in nearby clusters and 
associations. Although most of these do not reach the high 
accuracies realized for the Hyades, they hint at the future 
potential. 

Quantitatively, we have obtained radial velocities with 
standard errors of ~ 0.6 km s _1 for individual stars in the 
Hyades. The accuracies reached begin to make visible the 
convective and gravitational shifts expected in the spectra 
of F and G stars. For A stars and earlier types, where the 
convective shifts cannot yet be reliably predicted from the- 
ory, the spectra appear to be blueshifted by a few km s _1 
compared with the astromctrically determined motions 
and expected gravitational redshifts. This illustrates that 
astrometric radial velocities with uncertainties even well 
in excess of 1 km s _1 could be astrophysically interesting. 
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Such accuracies may also be sufficient to provide in- 
formation about the expansion of OB associations, as il- 
lustrated by the results for the Sco OB2 complex. Even 
with the modest precision of existing spectroscopic ve- 
locities, we see indications of expansion in the OB as- 
sociations Upper Centaurus Lupus and Lower Centaurus 
Crux (causing a bias in the astrometric radial velocities 
of 5-10 km s -1 ), while Upper Scorpius surprisingly shows 
no such indication. The limitations in the present under- 
standing of these associations come not from astrometry 
but mainly from spectroscopy and theory. 

From the same solution that gave astrometric ra- 
dial velocities, we get kinematically improved parallaxes. 
These can be used to study in greater detail the spatial 
structures and the Hertzsprung-Russell diagrams of both 
clusters and OB associations. In case of the Hyades, Upper 
Centaurus Lupus and Lower Centaurus Crux, the better- 
defined main sequences can also be taken as proof of the 
validity of the kinematic solution, and hence of the astro- 
metric radial velocities. 

Hipparcos parallax measurements reached typical ac- 
curacies of about 1.5 mas, while our improved parallaxes 
reach 0.3 mas for the Hyades. Space astrometry missions 
in the near future are expected to improve this by more 
than an order of magnitude to about 0.05 mas (Horner 
et al. 199§| ), with another order-of-magnitude gain by 
the future GAIA to 0.004 mas (Perryman et al. |200l| ). 
As detailed in Paper I, such accuracies will enable also 
other methods than the moving-cluster one for determin- 
ing radial velocities by purely geometric means. The future 
prospects for studying absolute radial velocities indepen- 
dent of spectroscopy look exciting indeed! 
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